Spectrum of Response to Platinum and PARP Inhibitors in Germline BRCA-Associated Pancreatic Cancer in the Clinical and Preclinical Setting.

Germline BRCA-associated pancreatic ductal adenocarcinoma (glBRCA PDAC) tumors are susceptible to platinum and PARP inhibition. The clinical outcomes of 125 patients with glBRCA PDAC were stratified based on the spectrum of response to platinum/PARP inhibition: (i) refractory [overall survival (OS) <6 months], (ii) durable response followed by acquired resistance (OS <36 months), and (iii) long-term responders (OS >36 months). Patient-derived xenografts (PDX) were generated from 25 patients with glBRCA PDAC at different clinical time points. Response to platinum/PARP inhibition in vivo and ex vivo culture (EVOC) correlated with clinical response. We deciphered the mechanisms of resistance in glBRCA PDAC and identified homologous recombination (HR) proficiency and secondary mutations restoring partial functionality as the most dominant resistant mechanism. Yet, a subset of HR-deficient (HRD) patients demonstrated clinical resistance. Their tumors displayed basal-like molecular subtype and were more aneuploid. Tumor mutational burden was high in HRD PDAC and significantly higher in tumors with secondary mutations. Anti-PD-1 attenuated tumor growth in a novel humanized glBRCA PDAC PDX model. This work demonstrates the utility of preclinical models, including EVOC, to predict the response of glBRCA PDAC to treatment, which has the potential to inform time-sensitive medical decisions. SIGNIFICANCE: glBRCA PDAC has a favorable response to platinum/PARP inhibition. However, most patients develop resistance. Additional treatment options for this unique subpopulation are needed. We generated model systems in PDXs and an ex vivo system (EVOC) that faithfully recapitulate these specific clinical scenarios as a platform to investigate the mechanisms of resistance for further drug development. This article is highlighted in the In This Issue feature, p. 1749.

hENT1 Expression Predicts Response to Gemcitabine and Nab-Paclitaxel in Advanced Pancreatic Ductal Adenocarcinoma.

PURPOSE: Modified FOLFIRINOX (mFFX) and gemcitabine/nab-paclitaxel (GnP) remain standard first-line options for patients with advanced pancreatic ductal adenocarcinoma (PDAC). Human equilibrative nucleoside transporter 1 (hENT1) was hypothesized to be a biomarker of gemcitabine in the adjuvant setting, with conflicting results. In this study, we explore hENT1 mRNA expression as a predictive biomarker in advanced PDAC. EXPERIMENTAL DESIGN: COMPASS was a prospective observational trial of patients with advanced PDAC. A biopsy was required prior to initiating chemotherapy, as determined by treating physician. Biopsies underwent laser capture microdissection prior to whole genome and RNA sequencing. The cut-off thresholds for hENT1 expression were determined using the maximal χ2 statistic. RESULTS: 253 patients were included in the analyses with a median follow-up of 32 months, with 138 patients receiving mFFX and 92 receiving GnP. In the intention to treat population, median overall survival (OS) was 10.0 months in hENT1high versus 7.9 months in hENT1low (P = 0.02). In patients receiving mFFX, there was no difference in overall response rate (ORR; 35% vs. 28%, P = 0.56) or median OS (10.6 vs. 10.5 months, P = 0.45). However, in patients treated with GnP, the ORR was significantly higher in hENT1high compared with hENT1low tumors (43% vs. 21%, P = 0.038). Median OS in this GnP-treated cohort was 10.6 months in hENT1high versus 6.7 months hENT1low (P < 0.001). In an interaction analysis, hENT1 was predictive of treatment response to GnP (interaction P = 0.002). CONCLUSIONS: In advanced PDAC, hENT1 mRNA expression predicts ORR and OS in patients receiving GnP.

Integrative analysis of KRAS wildtype metastatic pancreatic ductal adenocarcinoma reveals mutation and expression-based similarities to cholangiocarcinoma.

Oncogenic KRAS mutations are absent in approximately 10% of patients with metastatic pancreatic ductal adenocarcinoma (mPDAC) and may represent a subgroup of mPDAC with therapeutic options beyond standard-of-care cytotoxic chemotherapy. While distinct gene fusions have been implicated in KRAS wildtype mPDAC, information regarding other types of mutations remain limited, and gene expression patterns associated with KRAS wildtype mPDAC have not been reported. Here, we leverage sequencing data from the PanGen trial to perform comprehensive characterization of the molecular landscape of KRAS wildtype mPDAC and reveal increased frequency of chr1q amplification encompassing transcription factors PROX1 and NR5A2. By leveraging data from colorectal adenocarcinoma and cholangiocarcinoma samples, we highlight similarities between cholangiocarcinoma and KRAS wildtype mPDAC involving both mutation and expression-based signatures and validate these findings using an independent dataset. These data further establish KRAS wildtype mPDAC as a unique molecular entity, with therapeutic opportunities extending beyond gene fusion events.

Is Biannual Surveillance for Pancreatic Cancer Sufficient in Individuals With Genetic Syndromes or Familial Pancreatic Cancer?

BACKGROUND: Individuals with a family history of pancreatic adenocarcinoma (PC) or with a germline mutation in a PC susceptibility gene are at increased risk of developing PC. These high-risk individuals (HRIs) may benefit from PC surveillance. METHODS: A PC surveillance program was developed to evaluate the detection of premalignant lesions and early-stage PCs using biannual imaging and to determine whether locally advanced or metastatic PCs develop despite biannual surveillance. From January 2013 to April 2020, asymptomatic HRIs were enrolled and followed with alternating MRI and endoscopic ultrasound every 6 months. RESULTS: Of 75 HRIs, 43 (57.3%) had a germline mutation in a PC susceptibility gene and 32 (42.7%) had a familial pancreatic cancer (FPC) pedigree. Branch-duct intraductal papillary mucinous neoplasms (BD-IPMNs) were identified in 26 individuals (34.7%), but only 2 developed progressive lesions. One patient with Peutz-Jeghers syndrome (PJS) developed locally advanced PC arising from a BD-IPMN. Whole-genome sequencing of this patient's PC and of a second patient with PJS-associated PC from the same kindred revealed biallelic inactivation of STK11 in a KRAS-independent manner. A review of 3,853 patients from 2 PC registries identified an additional patient with PJS-associated PC. All 3 patients with PJS developed advanced PC consistent with the malignant transformation of an underlying BD-IPMN in <6 months. The other surveillance patient with a progressive lesion had FPC and underwent resection of a mixed-type IPMN that harbored polyclonal KRAS mutations. CONCLUSIONS: PC surveillance identifies a high prevalence of BD-IPMNs in HRIs. Patients with PJS with BD-IPMNs may be at risk for accelerated malignant transformation.

Intra-Tumoral CD8+ T-Cell Infiltration and PD-L1 Positivity in Homologous Recombination Deficient Pancreatic Ductal Adenocarcinoma.

The immune contexture of pancreatic ductal adenocarcinoma (PDAC) is generally immunosuppressive. A role for immune checkpoint inhibitors (ICIs) in PDAC has only been demonstrated for the rare and hypermutated mismatch repair (MMR) deficient (MMR-d) subtype. Homologous recombination repair (HR) deficient (HR-d) PDAC is more prevalent and may encompass up to 20% of PDAC. Its genomic instability may promote a T-cell mediated anti-tumor response with therapeutic sensitivity to ICIs. To investigate the immunogenicity of HR-d PDAC, we used multiplex immunohistochemistry (IHC) to compare the density and spatial distribution of CD8+ cytotoxic T-cells, FOXP3+ regulatory T-cells (Tregs), and CD68+ tumor-associated macrophages (TAMs) in HR-d versus HR/MMR-intact PDAC. We also evaluated the IHC positivity of programmed death-ligand 1 (PD-L1) across the subgroups. 192 tumors were evaluated and classified as HR/MMR-intact (n=166), HR-d (n=25) or MMR-d (n=1) based on germline testing and tumor molecular hallmarks. Intra-tumoral CD8+ T-cell infiltration was higher in HR-d versus HR/MMR-intact PDAC (p<0.0001), while CD8+ T-cell densities in the peri-tumoral and stromal regions were similar in both groups. HR-d PDAC also displayed increased intra-tumoral FOXP3+ Tregs (p=0.049) and had a higher CD8+:FOXP3+ ratio (p=0.023). CD68+ TAM expression was similar in HR-d and HR/MMR-intact PDAC. Finally, 6 of the 25 HR-d cases showed a PD-L1 Combined Positive Score of >=1, whereas none of the HR/MMR-intact cases met this threshold (p<0.00001). These results provide immunohistochemical evidence for intra-tumoral CD8+ T-cell enrichment and PD-L1 positivity in HR-d PDAC, suggesting that HR-d PDAC may be amenable to ICI treatment strategies.

Patient-derived tumor xenograft and organoid models established from resected pancreatic, duodenal and biliary cancers.

Patient-derived xenograft (PDX) and their xenograft-derived organoid (XDO) models that recapitulate the genotypic and phenotypic landscape of patient cancers could help to advance research and lead to improved clinical management. PDX models were established from 276 pancreato-duodenal and biliary cancer resections. Initial, passage 0 (P0) engraftment rates were 59% (118/199) for pancreatic, 86% (25/29) for duodenal, and 35% (17/48) for biliary ductal tumors. Pancreatic ductal adenocarcinoma (PDAC), had a P0 engraftment rate of 62% (105/169). KRAS mutant and wild-type PDAC models were molecularly profiled, and XDO models were generated to perform initial drug response evaluations. Subsets of PDAC PDX models showed global copy number variants and gene expression profiles that were retained with serial passaging, and they showed a spectrum of somatic mutations represented in patient tumors. PDAC XDO models were established, with a success rate of 71% (10/14). Pathway activation of KRAS-MAPK in PDXs was independent of KRAS mutational status. Four wild-type KRAS models were characterized by one with EGFR (L747-P753 del), two with BRAF alterations (N486_P490del or V600E), and one with triple negative KRAS/EGFR/BRAF. Model OCIP256, characterized by BRAF (N486-P490 del), had activated phospho-ERK. A combination treatment of a pan-RAF inhibitor (LY3009120) and a MEK inhibitor (trametinib) effectively suppressed phospho-ERK and inhibited growth of OCIP256 XDO and PDX models. PDAC/duodenal adenocarcinoma have high success rates forming PDX/organoid and retaining their phenotypic and genotypic features. These models may be effective tools to evaluate novel drug combination therapies.

Genomic Features and Classification of Homologous Recombination Deficient Pancreatic Ductal Adenocarcinoma.

BACKGROUND AND AIMS: Homologous recombination deficiency (HRD) in pancreatic ductal adenocarcinoma (PDAC), remains poorly defined beyond germline (g) alterations in BRCA1, BRCA2, and PALB2. METHODS: We interrogated whole genome sequencing (WGS) data on 391 patients, including 49 carriers of pathogenic variants (PVs) in gBRCA and PALB2. HRD classifiers were applied to the dataset and included (1) the genomic instability score (GIS) used by Myriad's MyChoice HRD assay; (2) substitution base signature 3 (SBS3); (3) HRDetect; and (4) structural variant (SV) burden. Clinical outcomes and responses to chemotherapy were correlated with HRD status. RESULTS: Biallelic tumor inactivation of gBRCA or PALB2 was evident in 43 of 49 germline carriers identifying HRD-PDAC. HRDetect (score ≥0.7) predicted gBRCA1/PALB2 deficiency with highest sensitivity (98%) and specificity (100%). HRD genomic tumor classifiers suggested that 7% to 10% of PDACs that do not harbor gBRCA/PALB2 have features of HRD. Of the somatic HRDetecthi cases, 69% were attributed to alterations in BRCA1/2, PALB2, RAD51C/D, and XRCC2, and a tandem duplicator phenotype. TP53 loss was more common in BRCA1- compared with BRCA2-associated HRD-PDAC. HRD status was not prognostic in resected PDAC; however in advanced disease the GIS (P = .02), SBS3 (P = .03), and HRDetect score (P = .005) were predictive of platinum response and superior survival. PVs in gATM (n = 6) or gCHEK2 (n = 2) did not result in HRD-PDAC by any of the classifiers. In 4 patients, BRCA2 reversion mutations associated with platinum resistance. CONCLUSIONS: Germline and parallel somatic profiling of PDAC outperforms germline testing alone in identifying HRD-PDAC. An additional 7% to 10% of patients without gBRCA/PALB2 mutations may benefit from DNA damage response agents.

Delving into Early-onset Pancreatic Ductal Adenocarcinoma: How Does Age Fit In?

PURPOSE: With the rising incidence of early-onset pancreatic cancer (EOPC), molecular characteristics that distinguish early-onset pancreatic ductal adenocarcinoma (PDAC) tumors from those arising at a later age are not well understood. EXPERIMENTAL DESIGN: We performed bioinformatic analysis of genomic and transcriptomic data generated from 269 advanced (metastatic or locally advanced) and 277 resectable PDAC tumor samples. Patient samples were stratified into EOPC (age of onset ≤55 years; n = 117), intermediate (age of onset 55-70 years; n = 264), and average (age of onset ≥70 years; n = 165) groups. Frequency of somatic mutations affecting genes commonly implicated in PDAC, as well as gene expression patterns, were compared between EOPC and all other groups. RESULTS: EOPC tumors showed significantly lower frequency of somatic single-nucleotide variant (SNV)/insertions/deletions (indel) in CDKN2A (P = 0.0017), and were more likely to achieve biallelic mutation of CDKN2A through homozygous copy loss as opposed to heterozygous copy loss coupled with a loss-of-function SNV/indel mutation, the latter of which was more common for tumors with later ages of onset (P = 1.5e-4). Transcription factor forkhead box protein C2 (FOXC2) was significantly upregulated in EOPC tumors (P = 0.032). Genes significantly correlated with FOXC2 in PDAC samples were enriched for gene sets related to epithelial-to-mesenchymal transition (EMT) and included VIM (P = 1.8e-8), CDH11 (P = 6.5e-5), and CDH2 (P = 2.4e-2). CONCLUSIONS: Our comprehensive analysis of sequencing data generated from a large cohort of PDAC patient samples highlights a distinctive pattern of biallelic CDKN2A mutation in EOPC tumors. Increased expression of FOXC2 in EOPC, with the correlation between FOXC2 and EMT pathways, represents novel molecular characteristics of EOPC.See related commentary by Lou, p. 8.

Subtype-Discordant Pancreatic Ductal Adenocarcinoma Tumors Show Intermediate Clinical and Molecular Characteristics.

PURPOSE: RNA-sequencing-based subtyping of pancreatic ductal adenocarcinoma (PDAC) has been reported by multiple research groups, each using different methodologies and patient cohorts. "Classical" and "basal-like" PDAC subtypes are associated with survival differences, with basal-like tumors associated with worse prognosis. We amalgamated various PDAC subtyping tools to evaluate the potential of such tools to be reliable in clinical practice. EXPERIMENTAL DESIGN: Sequencing data for 574 PDAC tumors was obtained from prospective trials and retrospective public databases. Six published PDAC subtyping strategies (Moffitt regression tools, clustering-based Moffitt, Collisson, Bailey, and Karasinska subtypes) were used on each sample, and results were tested for subtype call consistency and association with survival. RESULTS: Basal-like and classical subtype calls were concordant in 88% of patient samples, and survival outcomes were significantly different (P < 0.05) between prognostic subtypes. Twelve percent of tumors had subtype-discordant calls across the different methods, showing intermediate survival in univariate and multivariate survival analyses. Transcriptional profiles compatible with that of a hybrid subtype signature were observed for subtype-discordant tumors, in which classical and basal-like genes were concomitantly expressed. Subtype-discordant tumors showed intermediate molecular characteristics, including subtyping gene expression (P < 0.0001) and mutant KRAS allelic imbalance (P < 0.001). CONCLUSIONS: Nearly 1 in 6 patients with PDAC have tumors that fail to reliably fall into the classical or basal-like PDAC subtype categories, based on two regression tools aimed toward clinical practice. Rather, these patient tumors show intermediate prognostic and molecular traits. We propose close consideration of the non-binary nature of PDAC subtypes for future incorporation of subtyping into clinical practice.

A Preclinical Trial and Molecularly Annotated Patient Cohort Identify Predictive Biomarkers in Homologous Recombination-deficient Pancreatic Cancer.

PURPOSE: Pancreatic ductal adenocarcinoma (PDAC) arising in patients with a germline BRCA1 or BRCA2 (gBRCA) mutation may be sensitive to platinum and PARP inhibitors (PARPi). However, treatment stratification based on gBRCA mutational status alone is associated with heterogeneous responses. EXPERIMENTAL DESIGN: We performed a seven-arm preclinical trial consisting of 471 mice, representing 12 unique PDAC patient-derived xenografts, of which nine were gBRCA mutated. From 179 patients whose PDAC was whole-genome and transcriptome sequenced, we identified 21 cases with homologous recombination deficiency (HRD), and investigated prognostic biomarkers. RESULTS: We found that biallelic inactivation of BRCA1/BRCA2 is associated with genomic hallmarks of HRD and required for cisplatin and talazoparib (PARPi) sensitivity. However, HRD genomic hallmarks persisted in xenografts despite the emergence of therapy resistance, indicating the presence of a genomic scar. We identified tumor polyploidy and a low Ki67 index as predictors of poor cisplatin and talazoparib response. In patients with HRD PDAC, tumor polyploidy and a basal-like transcriptomic subtype were independent predictors of shorter survival. To facilitate clinical assignment of transcriptomic subtype, we developed a novel pragmatic two-marker assay (GATA6:KRT17). CONCLUSIONS: In summary, we propose a predictive and prognostic model of gBRCA-mutated PDAC on the basis of HRD genomic hallmarks, Ki67 index, tumor ploidy, and transcriptomic subtype.

Endogenous Retrovirus Transcript Levels Are Associated with Immunogenic Signatures in Multiple Metastatic Cancer Types.

Next-generation sequencing of solid tumors has revealed variable signatures of immunogenicity across tumors, but underlying molecular characteristics driving such variation are not fully understood. Although expression of endogenous retrovirus (ERV)-containing transcripts can provide a source of tumor-specific neoantigen in some cancer models, associations between ERV levels and immunogenicity across different types of metastatic cancer are not well established. We performed bioinformatics analysis of genomic, transcriptomic, and clinical data across an integrated cohort of 199 patients with metastatic breast, colorectal, and pancreatic ductal adenocarcinoma tumors. Within each cancer type, we identified a subgroup of viral mimicry tumors in which increased ERV levels were coupled with transcriptional signatures of autonomous antiviral response and immunogenicity. In addition, viral mimicry colorectal and pancreatic tumors showed increased expression of DNA demethylation gene TET2 Taken together, these data demonstrate the existence of an ERV-associated viral mimicry phenotype across three distinct metastatic cancer types, while indicating links between ERV abundance, epigenetic dysregulation, and immunogenicity.

GATA6 Expression Distinguishes Classical and Basal-like Subtypes in Advanced Pancreatic Cancer.

PURPOSE: To determine the impact of basal-like and classical subtypes in advanced pancreatic ductal adenocarcinoma (PDAC) and to explore GATA6 expression as a surrogate biomarker. EXPERIMENTAL DESIGN: Within the COMPASS trial, patients proceeding to chemotherapy for advanced PDAC undergo tumor biopsy for RNA-sequencing (RNA-seq). Overall response rate (ORR) and overall survival (OS) were stratified by subtypes and according to chemotherapy received. Correlation of GATA6 with the subtypes using gene expression profiling, in situ hybridization (ISH) was explored. RESULTS: Between December 2015 and May 2019, 195 patients (95%) had enough tissue for RNA-seq; 39 (20%) were classified as basal-like and 156 (80%) as classical. RECIST response data were available for 157 patients; 29 basal-like and 128 classical where the ORR was 10% versus 33%, respectively (P = 0.02). In patients with basal-like tumors treated with modified FOLFIRINOX (n = 22), the progression rate was 60% compared with 15% in classical PDAC (P = 0.0002). Median OS in the intention-to-treat population (n = 195) was 9.3 months for classical versus 5.9 months for basal-like PDAC (HR, 0.47; 95% confidence interval, 0.32-0.69; P = 0.0001). GATA6 expression by RNA-seq highly correlated with the classifier (P < 0.001) and ISH predicted the subtypes with sensitivity of 89% and specificity of 83%. In a multivariate analysis, GATA6 expression was prognostic (P = 0.02). In exploratory analyses, basal-like tumors, could be identified by keratin 5, were more hypoxic and enriched for a T-cell-inflamed gene expression signature. CONCLUSIONS: The basal-like subtype is chemoresistant and can be distinguished from classical PDAC by GATA6 expression.See related commentary by Collisson, p. 4715.

Author Correction: Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A Four-Chemokine Signature Is Associated with a T-cell-Inflamed Phenotype in Primary and Metastatic Pancreatic Cancer.

PURPOSE: The molecular drivers of antitumor immunity in pancreatic ductal adenocarcinoma (PDAC) are poorly understood, posing a major obstacle for the identification of patients potentially amenable for immune-checkpoint blockade or other novel strategies. Here, we explore the association of chemokine expression with effector T-cell infiltration in PDAC. EXPERIMENTAL DESIGN: Discovery cohorts comprised 113 primary resected PDAC and 107 PDAC liver metastases. Validation cohorts comprised 182 PDAC from The Cancer Genome Atlas and 92 PDACs from the Australian International Cancer Genome Consortium. We explored associations between immune cell counts by immunohistochemistry, chemokine expression, and transcriptional hallmarks of antitumor immunity by RNA sequencing (RNA-seq), and mutational burden by whole-genome sequencing. RESULTS: Among all known human chemokines, a coregulated set of four (CCL4, CCL5, CXCL9, and CXCL10) was strongly associated with CD8+ T-cell infiltration (P < 0.001). Expression of this "4-chemokine signature" positively correlated with transcriptional metrics of T-cell activation (ZAP70, ITK, and IL2RB), cytolytic activity (GZMA and PRF1), and immunosuppression (PDL1, PD1, CTLA4, TIM3, TIGIT, LAG3, FASLG, and IDO1). Furthermore, the 4-chemokine signature marked tumors with increased T-cell activation scores (MHC I presentation, T-cell/APC costimulation) and elevated expression of innate immune sensing pathways involved in T-cell priming (STING and NLRP3 inflammasome pathways, BATF3-driven dendritic cells). Importantly, expression of this 4-chemokine signature was consistently indicative of a T-cell-inflamed phenotype across primary PDAC and PDAC liver metastases. CONCLUSIONS: A conserved 4-chemokine signature marks resectable and metastatic PDAC tumors with an active antitumor phenotype. This could have implications for the appropriate selection of PDAC patients in immunotherapy trials.

Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution.

Pancreatic adenocarcinoma presents as a spectrum of a highly aggressive disease in patients. The basis of this disease heterogeneity has proved difficult to resolve due to poor tumor cellularity and extensive genomic instability. To address this, a dataset of whole genomes and transcriptomes was generated from purified epithelium of primary and metastatic tumors. Transcriptome analysis demonstrated that molecular subtypes are a product of a gene expression continuum driven by a mixture of intratumoral subpopulations, which was confirmed by single-cell analysis. Integrated whole-genome analysis uncovered that molecular subtypes are linked to specific copy number aberrations in genes such as mutant KRAS and GATA6. By mapping tumor genetic histories, tetraploidization emerged as a key mutational process behind these events. Taken together, these data support the premise that the constellation of genomic aberrations in the tumor gives rise to the molecular subtype, and that disease heterogeneity is due to ongoing genomic instability during progression.

Altered Gene Expression along the Glycolysis-Cholesterol Synthesis Axis Is Associated with Outcome in Pancreatic Cancer.

PURPOSE: Identification of clinically actionable molecular subtypes of pancreatic ductal adenocarcinoma (PDAC) is key to improving patient outcome. Intertumoral metabolic heterogeneity contributes to cancer survival and the balance between distinct metabolic pathways may influence PDAC outcome. We hypothesized that PDAC can be stratified into prognostic metabolic subgroups based on alterations in the expression of genes involved in glycolysis and cholesterol synthesis. EXPERIMENTAL DESIGN: We performed bioinformatics analysis of genomic, transcriptomic, and clinical data in an integrated cohort of 325 resectable and nonresectable PDAC. The resectable datasets included retrospective The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) cohorts. The nonresectable PDAC cohort studies included prospective COMPASS, PanGen, and BC Cancer Personalized OncoGenomics program (POG). RESULTS: On the basis of the median normalized expression of glycolytic and cholesterogenic genes, four subgroups were identified: quiescent, glycolytic, cholesterogenic, and mixed. Glycolytic tumors were associated with the shortest median survival in resectable (log-rank test P = 0.018) and metastatic settings (log-rank test P = 0.027). Patients with cholesterogenic tumors had the longest median survival. KRAS and MYC-amplified tumors had higher expression of glycolytic genes than tumors with normal or lost copies of the oncogenes (Wilcoxon rank sum test P = 0.015). Glycolytic tumors had the lowest expression of mitochondrial pyruvate carriers MPC1 and MPC2. Glycolytic and cholesterogenic gene expression correlated with the expression of prognostic PDAC subtype classifier genes. CONCLUSIONS: Metabolic classification specific to glycolytic and cholesterogenic pathways provides novel biological insight into previously established PDAC subtypes and may help develop personalized therapies targeting unique tumor metabolic profiles.See related commentary by Mehla and Singh, p. 6.

NRG1 Gene Fusions Are Recurrent, Clinically Actionable Gene Rearrangements in KRAS Wild-Type Pancreatic Ductal Adenocarcinoma.

PURPOSE: Gene fusions involving neuregulin 1 (NRG1) have been noted in multiple cancer types and have potential therapeutic implications. Although varying results have been reported in other cancer types, the efficacy of the HER-family kinase inhibitor afatinib in the treatment of NRG1 fusion-positive pancreatic ductal adenocarcinoma is not fully understood. EXPERIMENTAL DESIGN: Forty-seven patients with pancreatic ductal adenocarcinoma received comprehensive whole-genome and transcriptome sequencing and analysis. Two patients with gene fusions involving NRG1 received afatinib treatment, with response measured by pretreatment and posttreatment PET/CT imaging. RESULTS: Three of 47 (6%) patients with advanced pancreatic ductal adenocarcinoma were identified as KRAS wild type by whole-genome sequencing. All KRAS wild-type tumors were positive for gene fusions involving the ERBB3 ligand NRG1. Two of 3 patients with NRG1 fusion-positive tumors were treated with afatinib and demonstrated a significant and rapid response while on therapy. CONCLUSIONS: This work adds to a growing body of evidence that NRG1 gene fusions are recurrent, therapeutically actionable genomic events in pancreatic cancers. Based on the clinical outcomes described here, patients with KRAS wild-type tumors harboring NRG1 gene fusions may benefit from treatment with afatinib.See related commentary by Aguirre, p. 4589.

Meta-Analysis of 1,200 Transcriptomic Profiles Identifies a Prognostic Model for Pancreatic Ductal Adenocarcinoma.

PURPOSE: With a dismal 8% median 5-year overall survival, pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy. Only 10% to 20% of patients are eligible for surgery, and more than 50% of these patients will die within 1 year of surgery. Building a molecular predictor of early death would enable the selection of patients with PDAC who are at high risk. MATERIALS AND METHODS: We developed the Pancreatic Cancer Overall Survival Predictor (PCOSP), a prognostic model built from a unique set of 89 PDAC tumors in which gene expression was profiled using both microarray and sequencing platforms. We used a meta-analysis framework that was based on the binary gene pair method to create gene expression barcodes that were robust to biases arising from heterogeneous profiling platforms and batch effects. Leveraging the largest compendium of PDAC transcriptomic data sets to date, we show that PCOSP is a robust single-sample predictor of early death-1 year or less-after surgery in a subset of 823 samples with available transcriptomics and survival data. RESULTS: The PCOSP model was strongly and significantly prognostic, with a meta-estimate of the area under the receiver operating curve of 0.70 (P = 2.6E-22) and d-index (robust hazard ratio) of 1.9 (range, 1.6 to 2.3; ( = 1.4E-04) for binary and survival predictions, respectively. The prognostic value of PCOSP was independent of clinicopathologic parameters and molecular subtypes. Over-representation analysis of the PCOSP 2,619 gene pairs-1,070 unique genes-unveiled pathways associated with Hedgehog signaling, epithelial-mesenchymal transition, and extracellular matrix signaling. CONCLUSION: PCOSP could improve treatment decisions by identifying patients who will not benefit from standard surgery/chemotherapy but who may benefit from a more aggressive treatment approach or enrollment in a clinical trial.

Whole genomes define concordance of matched primary, xenograft, and organoid models of pancreas cancer.

Pancreatic ductal adenocarcinoma (PDAC) has the worst prognosis among solid malignancies and improved therapeutic strategies are needed to improve outcomes. Patient-derived xenografts (PDX) and patient-derived organoids (PDO) serve as promising tools to identify new drugs with therapeutic potential in PDAC. For these preclinical disease models to be effective, they should both recapitulate the molecular heterogeneity of PDAC and validate patient-specific therapeutic sensitivities. To date however, deep characterization of the molecular heterogeneity of PDAC PDX and PDO models and comparison with matched human tumour remains largely unaddressed at the whole genome level. We conducted a comprehensive assessment of the genetic landscape of 16 whole-genome pairs of tumours and matched PDX, from primary PDAC and liver metastasis, including a unique cohort of 5 'trios' of matched primary tumour, PDX, and PDO. We developed a pipeline to score concordance between PDAC models and their paired human tumours for genomic events, including mutations, structural variations, and copy number variations. Tumour-model comparisons of mutations displayed single-gene concordance across major PDAC driver genes, but relatively poor agreement across the greater mutational load. Genome-wide and chromosome-centric analysis of structural variation (SV) events highlights previously unrecognized concordance across chromosomes that demonstrate clustered SV events. We found that polyploidy presented a major challenge when assessing copy number changes; however, ploidy-corrected copy number states suggest good agreement between donor-model pairs. Collectively, our investigations highlight that while PDXs and PDOs may serve as tractable and transplantable systems for probing the molecular properties of PDAC, these models may best serve selective analyses across different levels of genomic complexity.